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Abstract 

Purpose Influenza A viruses, human coronaviruses 
(hCoV) and human bocavirus (hBoV) are emerging res¬ 
piratory viruses. This study investigated the association 
between influenza A viruses co-infection with hBoV and 
hCoV and severity and the sensitivity of a real-time poly¬ 
merase chain reaction (RT-PCR) assay for identification of 
15 coronaviruses. 

Methodology Published sequences for the 15 human 
coronaviruses were used to design a consensus PCR tar¬ 
geting the replicase open reading frame lb. A previously 
published PCR targeting the NS1 Gene of all known 
human bocavirus strains was also utilized. A series of 217 
samples from patients aged 37.7 (SD zb 30.4)] with sea¬ 
sonal influenza A viruses (SeasFluA) identified between 
06/2011 and 06/2012 in NW England were tested for 
hCoV and hBoV using RT-PCR. Association between co- 
infection and disease outcome was assessed using logistic 
regression. 

Results The limit of detection of hCoV RT-PCR assay 
was 2 copies/pl of human coronavirus RNA template, a 
sensitivity comparable to a previously published SYBR 
green assay for human coronaviruses. A total of 12 hCoV 
and 17 hBoV were identified in the 217 influenza A posi¬ 
tive samples. A higher proportion (61.5 %; 8/13) of 
SeasFluA/hBoV co-infections were identified in patients 
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that were admitted either to a general ward or the intensive 
care unit compared to 44.3 % (66/149) of single SeasFlu A 
virus infections (OR 2.5 95 % Cl 0.67-9.34, p = 0.17). In 
a stratified analysis, there was a trend towards higher asso¬ 
ciation between Flu A, hCoV and hBoV with increasing age 
(especially in patients aged 24-45 years and >65 year old). 
Conclusion Our hCoV RT-PCR protocol appeared to 
be of adequate analytical sensitivity for diagnosis. More 
and larger studies are needed to confirm the role of hCoV, 
hBoV in causing severe disease when they co-infect with 
influenza A viruses. 

Keywords Influenza A virus • Coronavirus • Bocavirus • 
SYBR Green /TaqMan RT-PCR assay • Co-infection • 
Hospitalization 

Introduction 

Influenza A viruses are a significant cause of acute lower 
respiratory tract infections (ALRI), with 13.0 % of all chil¬ 
dren with influenza infection developing ALRI-equivalent 
to 20 million (95 % Cl 13-32 million) ALRI cases per 
year [1]; and according to these authors, in 2008, between 
28,000 and 111,500 children younger than 5 years died of 
influenza-associated ALRI. Among potential co-pathogens, 
human coronaviruses (hCoV) and human bocavirus (hBoV) 
are important emerging respiratory viruses [2, 3]. 

Coronaviruses infect a wide variety of hosts (animals 
and humans); Alphacoronaviruses and Betacoronaviruses 
infect mammals, whereas Gammacoronaviruses infect 
birds [4]. Most of the human coronaviruses ( Alphacoro¬ 
naviruses 229E and NL63 and Betacoronaviruses OC43, 
HKU1) cause mild influenza like illnesses, however in 
2003, the severe acute respiratory syndrome (SARS) 
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coronavirus caused a very severe form of respiratory dis¬ 
ease in humans globally [5], resulting in over 900 deaths 
[6]. On 21 September 2012, a novel coronavirus was iden¬ 
tified in lower respiratory tract of an adult male Qatari 
national in London [7]. Three months earlier (13 June 
2012), Dutch researchers identified a similar virus in a 
Saudi Arabian national presenting with pneumonia [8-10]. 
As of 11 June, 2014, 699 laboratory-confirmed cases of 
human infection with MERS-CoV have been reported to 
WHO, including at least 209 deaths [6]. 

Human bocavirus, was first discovered in respira¬ 
tory specimens of patients with lower respiratory disease 
in 2005 by Allander et al. [11] using a random polymer¬ 
ase chain reaction (PCR)-cloning-sequencing protocol. 
Reviews by Schildgen et al. [3] and Zhenqiang et al. [12] 
showed that, since 2005, over 40 studies conducted glob¬ 
ally identified human bocavirus in children with acute res¬ 
piratory virus infections; and that most of these prevalence 
studies found hBoV occurring mainly together with other 
viruses. The significance of human bocaviruses in caus¬ 
ing serious clinical disease is still under debate, with some 
studies finding that bocavirus may cause significant ALRI 
[3, 12], and others finding they only cause mild illness [3, 
13, 14]. 

HBoV infection is reported to mainly occur as a co- 
infection with other respiratory viruses; some zoonotic 
hCoV have been associated with sever lower repiratory 
tract infection (LRTI). However these viruses are not, at 
present, routinely tested in clinical settings due to a lack 
of knowledge on the significance of these viruses in caus¬ 
ing severe disease. Some recent studies have indicated 
that respiratory virus co-infections increases disease 
severity [15, 16], however the pattern of co-infection 
between influenza A viruses and hCoV and hBoV and 
their association with severe disease has not been clearly 
elucidated. We designed primers and optimized a real¬ 
time polymerase chain reaction (RT-PCR) assay for iden¬ 
tification of 15 coronaviruses and used the RT-PCR tar¬ 
geting the NS1 region of the hBoV described by Qu et al. 
[17] to investigate the patterns of co-infections between 
influenza A virus with hCoV and hBoV and their associa¬ 
tion with severe clinical disease. 

Methodology 

Clinical samples and setting 

We tested for hCoV and hBoV in 217 samples randomly 
selected from 596 influenza A viruses identified from 3,804 
samples received between 24th June 2011 and 30th June 
2012 at the Manchester Microbiology Partnership Labora¬ 
tory (MMPL), a reference laboratory for the North West 


England serving a population of around 7 million peo¬ 
ple [18]. Influenza A viruses were identified using a well 
characterised in-house duplex real-time polymerase chain 
reaction assays and influenza A (H1N1) subtyped using 
the Public Health England (PHA) (Hl)v) assay [19]. The 
PCRs were run on the ABI7500 real-time PCR instrument 
(Applied Biosystems, Warrington, United Kingdom). Sam¬ 
ples were from patients aged between 1 and 98 years old 
(mean age 37.7, SD ± 30.4), seen as outpatients (at medi¬ 
cal centres, clinics or hospitals), or admitted to a hospital in 
North West England. 

Primers, templates and probes for hCoV, hBoV 

To design primers for hCoV, complete genomic sequences 
for 15 coronavirus (Table 1), were downloaded from the 
GenBank (and details of their GenBank accession num¬ 
bers available on request), aligned with BioEdit Sequence 
alignment Editor version 7.1.3.0 [20] and a consensus 
blasted using the BLAST program [GenBank, Bethesda 
MD, 20894 USA]; limiting the search to nucleotide col¬ 
lection (nt) database and organism Coronaviridae. One 
pair of primers targeting the region of the replicase open 
reading frame lb (ORFlb) non-structural protein 15 
(nspl5; XendoU/NendoU) uridylate specific endonu¬ 
clease, which is highly conserved among all coronavi¬ 
ruses [21, 22] was selected; forward primer: 5 7 -TGGGGA 
GTAATGAACCCGGTA-3' and reverse primer: S'-ACATGT 
A A A AG AGCTA ATA AC AC- 3' (Eurofins MWG Operon, 


Table 1 Coronaviruses used in the design of pan-coronavirus prim¬ 
ers (downloaded from GeneBank on 18th December, 2011) 



Coronavirus 

Number identified 
and downloaded 

1 

Human coronavirus HKU1 

20 

2 

Human coronavirus NL63 

2 

3 

Human coronavirus OC43 

4 

4 

Human coronavirus 229E 

1 

5 

Bovine coronavirus 

12 

6 

Porcine haemagglutinating virus 

2 

7 

Murine hepatitis virus 

21 

8 

Equine coronavirus 

3 

9 

Rat coronavirus 

1 

10 

Bat coronavirus 

31 

11 

Avian infectious bronchitis virus 

41 

12 

Canine coronavirus 

1 

13 

Sambar deer coronavirus 

1 

14 

Bat SARS coronavirus 

10 

15 

Duck coronavirus 

1 


Total 

151 
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Germany); probe 5 '-FAM-TT GT GTA A AG ATAGT GG AG 
GTT-MGB ~y (Applied Biosystems, Warrington, UK). 

The primer and template described by Qu et al. [17], 
GenBank accession no. NC_007455; forward primer 
S'-TAATGACTGCAGACAACGCCTAG-S', reverse primer 
S'-TGTCCCGCCCAAGATACACT-S' (Eurofins MWG 
Operon, Ebersberg Germany); probe 5' FAM-TTCCAC- 
CCAATCCTGGT-MGB-3' (Applied Biosystems, War¬ 
rington, UK), were used for identification of hBoV. 

Determination of analytical sensitivity and reproducibility 
of RT-PCR protocols 

To determine analytical sensitivity, synthetic DNA tem¬ 
plates were constructed using sequence nt position 
20062-20134 of Bovine coronavirus (accession number 
DQ811784.2 NCBI website): TGGGGAGTAATGAACCC 
GGTAATGTCGGTGGTAATGATGCTCTGGCAACCTC 
CACTATCTTTACACAAAGCCGTGTTATTAGCTCTTT- 
TACATGT and a sequence within the human bocavirus NS1 
gene: TAATGACTGCAGACAACGCCTAGTTGTTTGGT 
GGGAGGAGTGCTTAATGCACCAGGATTGGGTGG 
AACCTGCAAAGTGTATCTTGGGCGGGACAG. These 
were inserted into a pEX-A plasmid vector (Eurofins MWG 
Operon, Ebersberg Germany) and propagated in Escheri¬ 
chia coli. For use in the coronavirus RT-PCR, the plasmid 
was linearized using BamHl restriction enzyme (New Eng¬ 
land BioLabs, Ipswich USA) and RNA transcribed using 
the T7 High Yield RNA synthesis Kit (New England Bio- 
Labs, Ipswich USA). Purified RNA or DNA was quanti¬ 
fied using UV spectrophotometry and used to construct the 
standard curves used for quantitation in each of the PCR 
assays. To optimize the primer concentration, duplicate 
preparations of coronavirus 229E (obtained from Public 
Health England Type Culture Collection, Salisbury, UK), 
pEX-A RNA template, 20, 10, 7.50, 5 pM serial dilutions 
of forward and reverse primer and 10, 6.60, 3.3, 0.37 pM 
probe concentrations were set up in MicroAmp Fast Reac¬ 
tion PCR tubes and mounted on a 48 well plate. Similarly, 
for analytical sensitivity of the hBoV PCR, 1:4 serial dilu¬ 
tions (in duplicates) of the pEX-A plasmid (Eurofins MWG 
Operon, Ebersberg Germany) were used in standard curve 
experiments. 

PCR for identification of coronaviruses and bocavirus 
in samples 

The Power SYBR Green RNA-to-CT 1-Step Kit (Applied 
Biosystems) was used in a one-step RT-PCR for detection of 
coronavirus, the reaction mixture comprised: 10 pi of Power 
SYBR Green RT-PCR Master Mix (2x), 2.5 pi forward and 
reverse primer (20 pM), 1 pi of Array script RT Enzyme mix 
(125X), 2 pi of RNAse free water and 3 pi of RNA template 


in positive control well or sample in test wells. For identi¬ 
fication of bocavirus the Power SYBR Green PCR Master 
Mix (Applied Biosystems) was used: reaction mixture 10 pi 
of Power SYBR Green PCR Master Mix (2x), 2.5 pi for¬ 
ward and reverse primer (20 pM), 2 pi RNAse free water 
and 3 pi of template in positive control well or sample in test 
wells. RT-PCR experiments were conducted on the StepOne 
and StepOne Plus Real Time PCR machine (Applied Biosys¬ 
tems) with cycling parameters for coronavirus of: 48 °C for 
30 min, 95 °C for 10 min; followed by 40 cycles of dena¬ 
turing at 95 °C for 15 s and 57 °C for 1 min; and a melt¬ 
ing curve of 95 °C for 15 s, 57 °C for 15 s and 95 °C for 
15 s. The bocavirus PCR used 95 °C for 10 min, 40 cycles 
of 95 °C for 15 s and 60 °C for 1 min and a melting curve of 
95 °C for 15 s, 60 °C for 15 s and 95 °C for 15 s. The num¬ 
ber of viral copies was calculated using the absolute quantifi¬ 
cation method published by Brankatschk et al. [23]. 

Confirmation of positive samples using MGB probes 

Samples that were positive for coronavirus or bocavirus 
were confirmed with a TaqMan based RT-PCR. The probe 
concentrations were first optimized using duplicate prepa¬ 
rations of 10, 6.6, 5, 3.3, 0.3 pM and 0.12 probe concentra¬ 
tions. Each well contained 10 pi of TaqMan Fast Univer¬ 
sal PCR Master Mix (2X) (Applied Biosystems), 2.5 pi of 
probe, 2.5 pi forward and reverse primer, 2 pi of RNAse 
free water and 3 pd of RNA template. The coronavirus PCR 
also contained the 1 pi Multiscribe Reverse Transcriptase 
(Applied Biosystems) and 1 pi of RNase inhibitor. The 
cycling conditions were; for coronavirus experiments: 
50 °C for 2 min, 95 °C for 20 s, followed by 40 cycles of 
95 °C for 1 s and 57 °C for 20 s; whereas for bocavirus: 
95 °C for 20 s, followed by 40 cycles of 95 °C for 1 s and 
60 °C for 20 s. 

Polyacrylamide gel (PAGE) electrophoresis 

Amplicons for coronavirus and bocavirus PCR were visual¬ 
ized using precast 8 % TBE polyacrylamide gel 1 mm 15 
well (Applied Biosystems) run on an XCell SureLock™ 
Mini-Cell Electrophoresis chamber (Invitrogen Novex Mini 
cell), filled with 200 and 600 ml 1 x BTE Running Buffer 
in upper and lower chambers respectively. Ten (10 pi) 
of PCR product was mixed with 9 pi of Deionised water 
and 1 pi 5x DNA Loading Buffer (Bioline, London) and 
loaded into the wells and the positive control well loaded 
with HyperLadder V (Bioline). The loaded XCell Sure¬ 
Lock™ Mini-Cell chamber was run at 200 V, 18 mA for 
35 min, stained with GelRed Nucleic Acid Stain (Biotium, 
Hayward, UK) for 30 min, and visualised under UV light 
on an Alphaimager 2200 machine (Alpha Innotech, San 
Leandro, Calif). 
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21865 

21875 

21885 

21895 
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WYHYVHRGWH 

BWHSWNNBKK 

HRVYWKHBKN 

VKNYNHYRNH 

MBAGATTGGT 

CTCGGGYGTT 


21905 

21915 

21925 

21935 

21945 

21955 

Consensus 

CCATTAGACC 

CCTCATTACT 

TGGGCCATTA 

CAGCCACCAT 

TACTACGAGA 

CCGTTGGAKK 

F Primer 

-ACC 

CCTCATTACT 

TGGGCCATTA 

CA- 




21965 

21975 

21985 

21995 

22005 

22015 

Consensus 

BSDDAGWDRT 

SYDWHNYBDH 

RNDDHWWNHD 

WRNRRAWBYD 

SANVNDBNYD 

WWMNNWBHYD 

R Primer 








22025 

22035 

22045 

22055 

22065 

22075 

Consensus 

NYRDMMDWHN 

NVDWYNHVDB 

DHKKNRVRWD 

DVVBNNTBNV 

TRBYRNVVNT 

YYBDNYNHSV 


Fig. 1 A section of the consensus sequence (showing the region 
identical in all coronaviruses, derived from consensus sequences for 
15 coronaviruses) that was used to design pan-coronavirus primers. 
The sequence of region conserved in all coronaviruses is underlined , 

Statistical analysis 

Two by two tables were drawn to describe the patterns of 
co-infection between influenza A virus, human coronavirus 
and human bocavirus and differences in distribution of sin¬ 
gle and dual infections by age, sex, and season. Differences 
were assessed using the Pearson’s Chi square ( \ ) statistic, 
and for small samples, the Fisher’s exact test. Association 
between single and mixed virus infection and admission 
to a general ward (GW) or intensive care unit (ICU) were 
assessed using simple logistic regression models. All statis¬ 
tical analyses were performed with STATA software, ver¬ 
sion 11.0. (STATACorp, Texas 77845, USA). 

Ethics 

Ethical approval for the study was obtained from the 
National Health Service-Research Ethics Committee 
(NHS-REC), reference number ll/NW/0698), the Univer¬ 
sity of Manchester Research Ethics Office. Research and 
Development R&D approval was obtained from the Central 
Manchester Universities Hospitals NHS Foundation Trust 
(reference number R01835) (Fig. 1) 

Results 

Efficiency, analytical sensitivity and reproducibility of the 
hCoV and hBoV PCRs 

The optimum concentration for primers was 2.5 pi 
(20 pM) whereas the optimum concentration of the probe 


the positions here are in positive sense orientation, because of varia¬ 
tions in the sizes of coronavirus genomes, in bovine coronavirus, the 
amplicon region lies between nucleotide 20062 and 20156 


was 2.5 pi (6.6 pM). The optimum annealing temperature 
for the primers was 57 °C. At this temperature, the wells 
that contained coronavirus 229E and RNA template; gave 
the lowest threshold cycle (Ct value), (range 19-38), high¬ 
est fluorescence signal (range 0.6-19.2), and the PCR 
product melted at 76.2 °C, which is within the Tm for 
coronaviruses PCR products (75.0 and 80.0 °C). Analyti¬ 
cal sensitivity of the PCR assays was tested by seven 1:4 
serial dilutions (10 4 -10°) for the hCoV and eleven serial 
dilution (10 6 -10°) for the hBoV. The hCoV assay was able 
to detect down to 2.6 copies/pl, whereas the hBoV was 
able to detect down to 4.9 copies/pl. Under optimum con¬ 
ditions, the duplicate samples gave reproducible results 
and the experiments were repeated three times giving 
a total number of 21 and 33 replicates. Our results com¬ 
pare well with previously published protocols of hCoV 
and hBoV. For example a SYBR Green RT-PCR assay 
for identification of hCoV, published in 2007, Escutenaair 
et al. [24] reported that their assay was able to detect down 
to 10 RNA copies/pl, whereas Vijgen et al. [25] pan-cor- 
onavirus assay (published in 2008) was far less sensitive 
with a lower detection limit of 5.0 x 10 3 RNAcopies per 
microliter sample, Adach et al. [26] (1 and 10 genome 
copies). Further, hBoV assay by Neske et al. [27] was 
able to detect up to 3.8 x 10 s —3.8 X 10 3 copies/pl and Lu 
et al. [28] (10 copies). Our assays were therefore of good 
analytical sensitivity and could be used in diagnosis of 
hCoV. Polyacrylamide gel results indicated that the size 
of the PCR products obtained from the coronavirus 229E 
and bocavirus was similar to the coronavirus and bocavi¬ 
rus inserts into the pEX-A plasmid vector (a picture of the 
PAGE result for coronavirus is provided in Fig. 2). 
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Fig. 2 Polyacrylamide gel electrophoresis of hCoV PCR products. 
PCR products from two experiments, one where each well contained 
3 |xl of hCoV229E RNA (Expl) and another 5 |xl, to increase yield 
(Exp2) and products from RNA produced from the pEX-A vector. 
HypperLadder V marks are from bottom, 25, 50, 75, 100 bp up to 
500. The marks at 50 bp are due to primer dimmers which are known 
to form in SYBR Green PCRs 


Respiratory virus infections in samples that were positive 
for influenza A viruses 

Of the 217 influenza A virus samples included in this 
study, the majority 95.4 % (207) were seasonal influenza 
A viruses and only 10 (4.6 %) were pandemic influenza 
A(HlNl)pdm09 infections. This was because the pan¬ 
demic virus did not circulate in high numbers during the 
2011/2012 season. The World Health Organisation (WHO) 
had declared the end of the pandemic on 10th August 2010 
[29]. Human coronaviruses were found in 12 of the 217 
samples whereas human bocavirus in 17 of the 217 sam¬ 
ples. The Ct values for hCoV positive samples ranged 
from 21.95 to 35.0 (mean 31.3) corresponding to 122 and 
75 copies/pl respectively, whereas for hBoV it ranged from 
26.95 to 36.0 with (mean of 32.75) corresponding to 123 
and 92 copies/pl, indicating low viral loads. Previous stud¬ 
ies have also reported finding hBoV low and not high viral 
loads Zhenqiang et al. [12]. 

The same samples used for identification of hCoV 
and hBoV were also previously tested for other respira¬ 
tory virus infections: influenza B virus (Flu B), respira¬ 
tory syncytial virus (RSV), rhinoviruses (RV), adenovirus 
(AdV), human metapneumovirus (hMPV), and human 
parainfluenza virus types 1-3 (hPIVl-3). When the results 
of other respiratory viruses are included, 149 (72.0 %) of 


the 207 seasonal influenza A viruses were single infec¬ 
tions, and more than one virus was detected in 58 (28.0 %) 
cases. Only 1 of the 10 patients with pandemic influenza 
A(HlNl)pdm09 virus infection was also infected with 
human coronavirus and rhinovirus. Of the 34 co-infections 
between seasonal influenza A and other respiratory viruses; 
19 cases were infected with hMPV, 14 with RV, 2 with Flu 
B, 2 with RSV, and 2 with hPIV3 virus. 

Demographic and other characteristics of single and dual 
infections 

The age, sex, and season distribution of patients are given 
in (Table 2). The majority of seasonal influenza A single 
infections and hCoV and hBoV co-infections were in adults 
>19 years old. Specifically, 104/149 (69.8 %) of single sea¬ 
sonal influenza A virus infections were in adults >19 years 
compared to 8/11 (72.7 %) of SeasFluA/hCoV co-infec- 
tions and 9/13 (69.2 %) of SeasFluA/hBoV co-infections 
in adults >19 years old (and these differences were not 
statistically significant)-Table 2. However co-infections 
between Flu A and other respiratory viruses were signifi¬ 
cantly higher among children <5 years old (11.2 %-14/34; 
p = 0.03). Regarding seasonal distribution, most of both 
single and dual virus infections occurred during the win¬ 
ter season (with summer season having low virus activity). 
This is probably because we tested for hCoV and hBoV 
in samples that were positive for influenza A viruses and 
influenza viruses are known to circulate in high numbers 
mainly in winter seasons. 

A larger percentage (61.5 %; 8/13) of SeasFluA/hBoV 
co-infections were identified in patients that were admit¬ 
ted to either a general ward (GW) or the intensive care 
unit (ICU) compared to 44.3 % (66/149) of single Seas- 
Flu A virus infections (OR 2.5 95 % Cl 0.67-9.34, 
p = 0.17). In a stratified analysis, there was a trend towards 
higher association between Flu A, hCoV and hBoV with 
increasing age, but none of these were statistically sig¬ 
nificant (Table 2). In children <5 years, 50 % (1/2) of 
SeasFluA/hCoV and 66.7 % (2/3) of SeasFluA/hBoV co- 
infections were admitted to a GW or ICU compared to 
26.5 % (9/34) of single influenza A virus infections. Com¬ 
paratively, among patients aged 24-45 years old, 66.7 % 
(2/3) with SeasFluA/hCoV co-infection were admitted 
either to a GW/ICU compared to 56.3 % (18/32) of those 
with single seasonal influenza A virus infections (OR 
2.20, 95 % Cl 0.46-10.35, p = 0.32), and 100 % (3/3) of 
SeasFluA/hBoV co-infections were admitted to a GW/ICU. 
Similarly among those aged >65 years old, 75.0 % (3/4) of 
those with SeasFluA/hBoV co-infection were admitted to 
GW/ICU compared to 60.0 % (27/45) with single seasonal 
influenza A virus infections (OR 6.0, 95 % Cl 0.46-76.54, 
p = 0.17) and 100 % (1/1) of SeasFluA/hCoV co-infection 
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Table 2 Demographics and disease severity in single, dual and multiple seasonal influenza A virus infections 


Variable 

SeasFlu A only 

SeasFlu A + hCoV 

SeasFlu A + hBoV 

SeasFlu A + other RVIs 

Age (years )—n (%) 

<5 

34 (22.8) 

2(18.2) 

3 (23.1) 

14(41.2)* 

5-19 

11 (7.4) 

1 (9.1) 

1 (7.7) 

2(5.9) 

19-24 

6 (4.0) 

0 (0.0) 

0 (0.0) 

2(5.9) 

24 45 

32 (21.5) 

3 (27.3) 

3 (23.1) 

7 (20.6) 

45-65 

21 (14.1) 

4 (36,4 

2 (15.4) 

4(11.8) 

>65 

45 (30.2) 

1 (9.1) 

4 (30.8) 

5 (14.7) 

Total 

149(100) 

11 (100) 

13 (100) 

34(100) 

Male sex 

68 (49.3) 

7 (63.6) 

6 (46.2) 

17 (50.0) 

General ward/ICU —n (%) 

Age <5 

9/34 (26.5) 

1/2 (50.0) 

2/3 (66.7) 

3/14(21.43) 

>5 < 19 

0/0 (0.0) 

0/0 (0.0) 

0/0 (0.0) 

0/0 (0.0) 

>19 < 24 

1/6 (16.7) 

0/0 (0.0) 

0/0 (0.0) 

0/2 (0.0) 

>24 < 45 

18/32 (56.3) 

2/3 (66.7) 

3/3 (100) 

4/7 (57.1) 

>45 < 65 

11/21 (52.4) 

0/4 (0.0) 

0/2 (0.0) 

2/4 (50.0) 

>65 

27/45 (60.0) 

1/1 (100) 

3/4 (75.0) 

3/5 (60.0) 

Total 

66/149 (44.3) 

4/11 (36.6) 

8/13 (61.5) 

12/34 (35.3) 

OR (95 % Cl; p value) a 

Baseline 

0.60 (0.15-2.40; 0.47) 

2.5 (0.67-9.34; 0.17) 

0.87 (0.37- 2.10; 0.76) 


SeasFlu A seasonal influenza A virus, hCoV human coronavirus, hBoV human bocavirus, other RVIs co-infection between seasonal influenza A 
and other respiratory viruses (viruses involved; hBoV = 4, rhinovirus = 14, influenza B = 2, respiratory syncytial virus = 2, adenovirus = 4, 
human metapneumovirus =19 and parainfluenza virus type 3 = 1) 

The statistic was significant at p < 0.05 
a Odds ratios were adjusted for age group and sex 


was admitted to a GW/ICU (Table 2). The lack of statisti¬ 
cal significance was probably because most of the cells in 
the >5-19, >19-24 and >45-65 years age groups had zero 
observations; when these age groups were eliminated from 
the analysis, SeasFluA/hBoV co-infections was associated 
with a significant increase in risk of admission to GW/ICU 
(OR: 5.83, 95% Cl: 1.03-32.93,/? = 0.05). 

Discussion and conclusion 

The importance of developing a coronavirus RT-PCR assay 
that is able to identify both Alphacoronaviruses and Betac- 
oronaviruses which attack animals and humans, and Gam- 
mac oronaviruses which attack birds, is that such assays 
could identify zoonotic coronaviruses associated with more 
severe disease in humans [2, 5, 30, 31]. Therefore the clini¬ 
cal usability of any PCR assay for the identification of all 
or most coronaviruses depends on its ability to amplify a 
region which is common to all coronaviruses (both human 
and human adapted animal viruses), and its sensitivity. 
The primers for coronavirus we designed target the repli- 
case ORFlb nspl5 (XendoU/NendoU) uridylate specific 
endonuclease gene, which is highly conserved among all 
coronaviruses [21, 22] and we confirmed our results with 
a TaqMan probe. The quality of primers depends, among 


others, on their ability not to give nonspecific amplifi¬ 
cations. Our primers for identification of coronaviruses 
did not give nonspecific amplifications as indicated by 
the melting curves. In addition, the hCoV standard curve 
experiments indicated that the hCoV RT-PCR assay was 
able to pick up to 2 copies/pl of human coronavirus RNA 
template, our assays results compare equally with those of 
a previously published SYBR green assay for coronavi¬ 
ruses [24-26]. We are therefore confident that it could be 
applied for diagnosis of coronaviruses in clinical settings 
subject to further experiments in statistically significant 
number of samples. 

Some studies have reported that co-infections are pre¬ 
dominantly in younger and older patients [32-34]. In this 
study, single seasonal influenza A viruses were predomi¬ 
nantly in adults >18 years old in agreement with our ear¬ 
lier published study [35]. Similarly a higher proportion 
of hCoV and hBoV co-infections were in adults 18 to 
65 years old. This suggests high prevalence of hCoVs and 
hBoV among older patients. Comparatively the majority of 
other respiratory virus co-infections (RSV, RV, hMPV and 
hPIVl-3) were in children <5 years old (Table 2). Regard¬ 
ing seasonality of the viruses, almost all the respiratory 
virus single and dual infections in this study were identi¬ 
fied in the winter season. This could be because we tested 
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for hCoV and hBoV in samples that were positive for influ¬ 
enza A viruses and influenza viruses are known to circu¬ 
late mainly during winter season. Our findings are in agree¬ 
ment with other studies [16, 36], which also reported dual 
or multiple respiratory virus infections mainly occur during 
winter season. 

This study did not find a significant association between 
seasonal influenza A viruses co-infection with hCoV and 
hBoV and severe disease. Given the high co-infection rates 
for bocavirus reported by other studies [3, 12], probably 
more bocavirus infections could have been identified if all 
the 3,804 samples received at the MMPL between June 
2011 and June 2012 were tested. A possibility of finding 
more coronaviruses could also not be ruled out. A review 
by Schildgen et al. [3] suggested that human bocavirus may 
be a bystander and not a real pathogen; one of the evidence 
to this hypothesis being frequently identified co-infecting 
with other respiratory viruses, and also it is mostly iden¬ 
tified in low copy numbers. However, the results of this 
study indicate a trend towards the hypothesis that hBoV co- 
infection might promote hospitalization especially among 
patients >65 years old and >24 <45 years old. Further, a 
number of other factors; age, comorbidities, immune status, 
and bacterial co-infection, affect disease outcome [37-39]. 
However in this study we did not control for these factors 
as the datum did not have these variables. This may have 
affected our outcomes and this should be born in mind 
when interpreting our results. 

In conclusion, our RT-PCR protocol for the identification 
of coronaviruses is sensitive, specific and reliable and could 
be used for diagnosis of coronaviruses. More and larger stud¬ 
ies are needed to confirm the role of hBoV in causing severe 
disease when they co-infect with influenza A viruses. 
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